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We review the perturbative evolution of the polarized structure functions 
gi and their associated parton distribution functions, with particular em- 
phasis on the anomalous coupling of the first moment of the polarized 
gluon distribution. We also describe the small x behaviour of polarized 
parton distributions, contrasting it with that of the unpolarized distribu- 
tions. We then explain how this theoretical analysis affects the extraction 
of the singlet axial charge from experimental data on gi , and show that it 
may be possible to use such data to infer the existence of polarized gluons 
in the nucleon. 



The publication of the EMC results^ for the polarized proton structure func- 
tion has been directly responsible for a renewed interest in polarized deep 
inelastic scattering among the theoretical community. In particular the im- 
plication that in the naive parton model the total helicity carried by quarks 
and antiquarks in the proton was consistent with zero led to a rexamination of 
the role played by the axial anomaly and polarized gluons in the perturbative 
evolution of the first moment of polarized structure functions.^"'' Since then, 
more precise measurements^" of both g\ and gf over a wider range of both 
X and have been made, and from them seemingly very precise values of 
the first moments deduced. ^^'^^ These first moments are generally obtained 
by extrapolating the experimental data to a common scale (which is done in 
practice by assuming that the asymmetries measured in the experiments are 
independent), and further by extrapolating from the measured region to 
small X (using Regge behaviour, and assuming that the small x contribution is 
then Q"^ independent). However if the polarized gluon distribution were large, 
both of these approximations could turn out to be very poor,^** because of 
the anomalously large coupling of polarized gluons to the first moment of gi . 
It thus becomes necessary to examine in detail the theoretical errors implicit 
in our present ignorance of the size of the polarized gluon distribution, and 
conversely whether the x and dependence of existing or future structure 
function data may be used to infer the existence of a large gluonic contribution 
to the nucleon spin. 
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1 Polarized Partons 



Wc begin by reviewing the relation between polarized structure functions 
and parton densities in the parton model, and its relation to the operator 
product expansion and renormalization group approach. In the parton model^^ 
the polarized structure function gi is decomposed in terms of polarized quark 
and gluon distributions Aqi and Ag according to 



9i{x,Q^) /'^ |cNs(f ,a«(t))A(7Ns(2/,t) 

^ J X y ^ 



(1.1) 

+Cs(|, a,(f))Ags(2/, t) + C,(|, a,{t))Ag{y, t)} + 0{l/Q') , 

where t = \n^, the various coefficient functions C{x,as) are directly related 

to hard cross-sections calculable in perturbative QCD, and A^ns and Aqs are 
respectively the nonsinglet and singlet quark distributions: 



Aq^six,t) = ^^^{Aqi{x,t) + Aqi{x,t)), 



i=l 



(1.2) 



Aqs{x, t) = Y^{Aqi{x, t) + Aqi{x, t)), 



i=l 



where nf is the number of active flavours, each with electric charge e^, and 
(e^) = '^ef/uf. Although they are themselves intrinsically nonperturbative, 
the perturbative part of the x and t dependence of the polarized quark and 
gluon distributions is given by Altarelli-Parisi equations:^^ the nonsinglet quark 
evolves independently as 



j^Aq^s{x,t) = ^ £^Pf/(|,a,(t))AgNs(2/,i), (1-3) 



while the singlet quark and the gluon mix according to 

lAqs{x,t) = ^j'^ [p^^i^,a,{t))Aqs{y,t)+P,,i^,a,it))Agiy,t) 

±Ag{x,t) = ^ £^ [p,,(|,a,(t))Ags(y,t)+P,,(f,a,(0)A5(2/,t) 

(1-4) 

The splitting functions P{x, ag) are again computable perturbatively in terms 
of hard cross-sections. In the naive parton model Cns = C's = 5{l — x), Cg = 0, 
so the gluons decouple. In LO perturbation theory the gluons couple through 
the singlet evolution equations (1.4), while at NLO they also couple directly. 
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Heavy quark contributions are generated radiatively as thresholds are crossed. 
All this works in just the same way as in the unpolarized case. 

Taking the Mellin transform of eq.(l.l) gives the leading twist component 
of the operator product expansion of the moments of gi : 

= ^ [C^AqW + C%Aq% + Cf,Ag^] + 0{l/Q'). 

Here CAr(as) = /^^ dxx^^^C{x,as) are the Wilson coefficients and AqN{t) = 

/p dxx^-'^Aq{x,t) may, for some values of N, be related to forward matrix el- 
ements of local operators. In the unpolarized case moments of the distributions 
qi + qi and g correspond to matrix elements of local twist two operators for 
N = 2,4:, . . ., while for qi — qi suitable local operators exist for A'' = 1, 3, . . .: all 
the other moments are well defined, but can only be obtained by analytic con- 
tinuation in N (a good example^'' being the Gottried sum qf^). Here however 
the situation is reversed: 



^flPfl • • • PvN-l 

AgN{t) 



= {p,s\ol::^u-,^,jp,^)t, 



TV = 1,3,5,... 

TV = 3,5,7,... (1.6) 

Ar = 3,5,7,..., 



where \p, s) is some hadronic state carrying momentum p, with polarization 
vector Sfi, and the twist two local operators 0% are purely gluonic. Although 
there exist local operators for even N = 2,4, . . ., these have opposite charge 
conjugation and thus correspond to moments of the valence distributions qi — qi- 
Again all other moments are well defined, but can only be obtained by analytic 
continuation. They are necessarily gauge invariant since the matrix elements 
(1.6) are gauge invariant, but will in general be scheme dependent. The case 
A = 1 in the singlet channel is peculiar, in that there is just one local operator, 
the axial singlet current: the identification of matrix elements of this operator 
with the first moments Aqf and Agi is then rather subtle and will be discussed 
in the next section. 

Taking moments of the Altarelli-Parisi equations (1.3), (1.4) yields the 
renormalization group equations for the matrix elements (1.6): 

^A4^(t)=7^-^^(a.(t))A4^(t) 

l[Aq%{t)\^fj%^'^as{t)) 7^n«s(t))WAd(t)\ ^^-^^ 
dt { AgM it) J { (a, (t)) {a, (t)) J { Ag^ (t) ) ' 

where ^n{ois) = /g^ dxx^~^ P{x,as) are the anomalous dimensions of the 
various local operators. 
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2 First Moments 

Due to interest in the Bjorken^^ and Ellis-Jaffe^" sum rules, much of the 
hterature on gi{x,Q'^) has focussed on its first moment ri{Q^). The proper 
interpretation^"^ of the singlet first moments is complicated by the presence 
in this channel of the axial anomaly.^^ 

Consider firstly the renormalization group equations (1.7) with A^" = 1. 
While some elements of the matrix of anomalous dimensions are nonvanishing 
at LO (one loop), this dependence turns out to be trivial, since the eigenvectors 
of the evolution 

Aqf^, Aqf and ao = Aqf - Uf^Ag^, (2.1) 

only evolve at NLO (two loops). ^ It follows that there must exist factorization 
schemes in which both Aqf^ and Aqf do not evolve at all, since multiplicative 
renormalizations which only begin at NLO can always be removed by a change 
of scheme. Thus in such schemes 

|Agf=o, iAqf=0 and ^ao = 7.«o, (2.2) 

where 7s = —71/(0:^/277^) H has been calculated at two^^ and three^^ loops. 

In the nonsinglet sector such schemes are essential in order to make the 
usual identification (1.6) of Agf ^ with forward matrix elements of (partially) 
conserved nonsinglet axial currents, i.e. with the nonsinglet axial charges. 
These may then be determined (assuming exact flavor symmetry) from weak 
decays of hyperons: below the charm threshold 

^qf^ = + iag, (2.3) 

where Sf^gA = {p,s\{J^)3\p,s), etc. Above heavy quark thresholds nonsinglet 
charges ais, a24, etc. must also be added. The nonsinglet combinations gA = 
Au — Ad and as = Au + Ad — 2As, etc. of the quark contributions to the 
spin of the hadron are then both well defined and scale independent (above 
threshold). The Bjorken sum rule^^ follows immediately from eqn.(l.-5). 

Similarly in the singlet sector the conservation of Aqf makes it the natural 
candidate^" for the singlet quark contribution Au + Ad + As + ■ ■ ■ to the 
hadron spin. Individual quark contributions may then be disentangled, and 
in particular the Zweig rule ~ As ~ Ac ~ • • • acquires a scale independent 
meaning. In fact although there is no local gauge invariant conserved current 
in the axial singlet channel, Aqf may still be formally identified with singlet 
quark helicity.^'* The other evolution eigenvector uq may then be identified with 
forward matrix elements {p, s\j^\p, s)t = s^ao(t) of the singlet axial current 
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j^, whose conservation is violated at one loop by the (purely gluonic) axial 
anomaly, 

dfijl, — nf^e^j^i,paUG^i,Gp„, (2.4) 

which is thus directly responsible for the NLO evolution (2.2) of ao- Since the 
anomaly (2.4) is unaffected by higher order perturbativc corrections,^^ it is 
possible to find schemes in which both the decomposition (2.1) of and its 
evolution (2.2) hold to all orders in perturbation theory:^ in such schemes 

= ll\as) =-^s{as) + (i{ots), (2.5) 

where P{as) = d\naa/dt. 

Now consider the Wilson coefficients in the operator product expansion 
(1.5). When these are calculated in a factorization scheme in which first mo- 
ments evolve according to (2.2), their first moments are 

Cr = l-^ + ---, Cf = l-^ + ---, Cf = -nf^ + ---, (2.6) 

at NLO (in fact Cf^ is known to^e 0{al), Cf to^^ 0{al)). When combined 
with (2.1), this implies that when N = 1 there are indeed only two terms on 
the right hand side of (1.5): 

ri(Q') = ^ [Cf^{a,)Aqf^{t) + CfiaMt)] , (2-7) 

in accordance with the fact that there is only one local gauge invariant singlet 
operator with twist two and spin one, the axial singlet current j^, which must 
then be multiplicatively renormalized. Although (1.5), (2. 6) and (2.1) only 
imply (2.7) at NLO, in schemes in which the Adler-Bardeen theorem is satisfied 
it must be true to all orders in perturbation theory: in such schemes 

Cf(a,)=-n/fjCf(a,). (2.8) 

In practice all of these results may be obtained by regularizing infrared 
collinear divergences by putting external particles off-shell.^* Alternatively the 
infrared divergences may regulated by giving the quarks a mass;^ this has the 
disadvantage that chiral symmetry is broken, so a finite renormalization must 
be performed to ensure that the quark helicities A^/J^^ and Aqf remain scale 
independent. More seriously, if dimensional regularization alone is used, with 
masslcss qiiarks and all external particles on-shell, Cf vanishes, so Aqf is iden- 
tified with the axial singlet charge and thus can no longer be related to quark 
helicity. Such factorization schemes are inappropriate because soft contribu- 
tions are included in the coefficient function, rather than being factorized into 
the parton densities. 
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Finally wc consider the implications of this analysis. Combining (2.1) 
and (2.2), we see that Agx{t) ~ f as t — » cxd, so asymptotically the gluonic 
contribution to ao{t) does not decouple, but may be as large as the quark 
contribution.^ Thus even if Aqf ^ ag as suggested by the Zweig rule, o,o{t) 
may still be very different from as, and the Ellis- Jaffe sum rule^° may fail. 

To actually explain why experimentally^ ao(t) is small in the perturbative 
region we would nevertheless require some nonpertm-bative mechanism. Vari- 
ous attempts at natural expanantions exist, in which either Aqf is suppressed 
by instantons,^^ and thus the Zweig rule is strongly violated, or else the Zweig 
rule holds but ao{t) is suppressed, either by strong evolution at low scales, or 
due to the smallness of the first derivative of the topological susceptibility.^" 
The latter explanations are both target independent, so might be tested by 
measuring gi for other targets such as the photon. Meanwhile, it would be 
useful to have some way of determining Aq^ or A^f independently. 

Since the gluonic contribution to the first moment of gi is effectively LO, it 
might be reasonable to expect that if Ag were indeed large the scale dependence 
of gi{x, Q"^) might also be anomalously large in certain regions of x}'^ If this 
were so it might be possible conversely to determine A^ by studying scaling 
violations of gi. Wc will return to this idea in section 4, but first will discuss 
another way of arriving at the same conclusion. 



3 Small X Evolution 

The behaviour of structure functions at small x is governed both by the 
qualitative form of the non-perturbative input at some scale to = InQo/A^, 
given in principle by Regge theory, and the structure of the perturbative evo- 
lution to > Qq, which can then often overwhelm it. Indeed the form 
and evolution of structure functions at small x depends on the outcome of a 
competition between Regge behaviour, LO evolution and higher order pertur- 
bative evolution. We will begin by describing how this works for unpolarized 
distributions, and then consider the polarized distributions. 

According to Regge theory the behaviour as a; —> of the unpolarized 
singlet distributions is controlled by the pomeron intercept, so xqs{x,to) ~ 
xg{x,to) ^ x~'^, with A ~ 0.08. However this essentially flat (or 'soft') be- 
haviour is substantially modified by perturbative evolution, determined by the 
leading singularity which for singlet anomalous dimensions is at = 1 . At one 
loop, 7^ ~ 7^ ~ as/ {N—1) as N 1, while 7^ and 7^ are both regular. It 
follows that in the double limit a; — > and f — > 00 the gluon distribution grows 
in a precisely determined way, faster than any power of In 1/a; but slower than 
any power oi x:^^ 

xg{x, t) ~ AAcr-i/^ exp(27o- - SQ, (3.1) 
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where if ^ = In = In ^ In a = \/^, and 7 and S arc (known) 

numerical constants. This growth in xg drives a similar 'double scaling' be- 
haviour in xqs , which has recently been confirmed by measurements of F2 at 
small X and large t at HERA.^^ 

Two loop corrections to double scaling are small^^ essentially because 
the singularity in the two loop singlet anomalous dimensions is of the form 
a1/{N — 1), and thus no stronger than that at one loop. Naively one expects 
that at / loops 7!^^ ~ ~ because at each order there is both an 

extra mass singularity and an extra collinear singularity,^^ which would wreck 
double scaling by inducing a strong power-like growth. However many of the 
singularities cancel, '^^''^^ so that the true behaviour at / loops is 7^ ^ {as/{N— 
1))'. The (scheme independent) coefficients of these remaining singularities can 
be calculated'^^ and in fact turn out to be very small; indeed the cocfhcients 
at 2, 3 and 5 loops actually vanish. This means that although the higher 
singularities may eventually produce a power-like growth in the Regge limit 
like where Ag = 12ln2as/TT is the radius of convergence of the sum of 

singularities, no indication of such behaviour has yet been seen at HERA.^^'^^ 

The behaviour of nonsinglet unpolarized parton distributions at small x 
is rather different, however. To begin with, the leading singularity in 7^'^^^ is 
not at = 1 but at A^ = 0: all singularities at A^ = 1 cancel. According to 
Regge behaviour as a; — > qNs{x,to) ^ a;""^, with A ~ 0.5, which now counts 
as a 'hard' boundary condition. So instead of the double scaling behaviour 
(3.1), we now have 

QNsix, t) ~ A'exp(A^ + (f/A - 5)Q, (3.2) 

for p = \/ljC ^5 7/ A, 7 and 5 constants. It follows immediately that at small 
X qns is considerably smaller qs, and thus — will be much more difficult 
to measure than alone. The two loop correction is now rather larger since 
to 0{al) jf^'^^ behave as a1/N^ as N ^ 0, and indeed at I loops there is now 
no cancellation of double logarithmic singularities,^^ so 7^' ~ a'/A"^'"^. 
Again the coefficients of the singularities can be computed,'"' and turn out 
to be large. ^ It follows that even if qNsi^-U)) had been soft, asymptotically 
QNs{x,t) will eventually grow as x~^^^ , where^^ Ans = yoOs/STT ~ |. As 
it is, this behaviour should set in very quickly, and indeed it seems to be in 
agreement with existing NMC and CCFR data. 

* Although there is now no factorization theorem to gaurantee that the resulting leading 
singularity anomalous dimension is scheme independent, the results at two loops (both un- 
polarized and polarized''^) are correct,''^ and the behaviour of the MS coefficient functions at 
0{as) and*^ 0(al) is not so singular as to introduce scheme dependence into the anomalous 
dimension (at least at NNLO). It thus seems that the imphcations of ref.39 deserve to be taken 
seriously. 
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Polarized distributions, both nonsinglet and singlet, have a soft Rcgge 
behaviour:''^'"'^ Agjvs(a;,to) ~ ^Qs{x,to) ~ x~^, with —0.5 ^ A ^ while 
again all perturbative singularities a.t N = 1 cancel. Indeed at LO all anoma- 
lous dimensions behave as ag/N, and thus one might expect a double scaling 
growth for both nonsinglet and singlet distributions: ^"^'^^"^^ 



where 7ns > 7±) ^ns and 6± are all (known) constants, and v±{x,t) are eigen- 
vectors of the evolution: v± = {Aq^,Ag^), Aq^ = —c±Ag^, with c± both 
positive constants. Thus rather than the gluon driving the singlet quark, as 
happens in the unpolarized case, here both gluon and singlet quark grow to- 
gether, but with opposite signs. It follows that in general while A^ns and 
Ag become large and positive as a; ^ and t ^ oo, Aqs becomes large and 
negative, and gi{x, Q"^) can then exhibit strong fluctuations due to interference 
between the various contributions.^^ Furthermore g\ and can behave rather 
differently, since AgNS grows on the same footing as Ags- 

However, these results must be interpreted with great care since at small 
enough x higher order singularities can quickly become important. All the two 
loop anomalous dimensions behave as a^/N^, and although the two loop cor- 
rections are still numerically small in the small x region accessible currently,^^ 
and in particular leave the singlet eigenvectors unchanged at NLO, at yet 
smaller x they soon become comparable to the LO term. Just as in the non- 
singlet unpolarized case, there is no evidence that the double logarithmic sin- 
gularities cancel: indeed the leading singularities in the nonsinglet polarized 
channel have already been computed,'''' and the resulting power behaviour 
of Aqi^s{x,t) turns out to be even stronger^° than that of qj<!s{x,t). The 
same seems to be true in the polarized singlet channel, though the matrix of 
anomalous dimensions is not yet known. 

To summarise, as 1/x and increase grows rapidly due to 

perturbative evolution, and may go negative if Ag is large enough. The precise 
details of the behaviour at very small x have yet to be calculated, but will 
probably involve rather complicated fluctuations in the {x, Q^) plane. However 
what is already clear is that the small x contribution to the first moment of gi 
may have a strong dependence. Since the scale dependence of the complete 
first moment (2.7) is perturbatively rather weak, this would necessarily imply 
a compensating dependence at larger x driven by Ag, which may be visible 
in existing data. 

4 Polarized Gluons? 

Fixed target experiments with a fixed beam energy measure the polariza- 
tion asymmetries, and thus indirectly gi{x,Q'^), along a curve Ql {x). For 



Aqy!s{x,t) ~ A/kso- 

V±{x,t) ~ M±(T~ 



1/2 



^ exp(27Nso- - (5nsC) 
exp(27±cr - (5±C), 



-1/2 



(3.3) 
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dashed: Q^=Q^(x) (SMC) 
solid: Q^=Q^(x) (E143) 
dotdashcd: q''=1 GeV* 
dotted: q^=10 GeV^ 




(i) b 




Figure 1: The SMC*'" (crosses) and EUS^O'ii (diamonds) data for (a) and (b) gf plotted 
against x. The curves correspond to (i) maximal gluon and (ii) minimal gluon flts.^* 



the SMC experiments*^'^ this curve goes from (0.7, 50GeV^) to (0.003, IGeV^), 
while for the E143 experiments the beam energy is fower, so the is fower 
for aU X, the curve reaching from (0.7,9GeV^) to (0.03, IGeV^). Combining 
both experiments thus gives us two values of for each value of x, making it 
possible to search for purely evolutionary effects. 

Several fits to the available data have been made, by evolving a standard 
parameterization of A(jNs(ic, to); ^Qsi^ito) and Ag{x,to) from the starting 
scale (usually taken to be Qq = IGeV^) up to the data. Until recently the 
evolution was performed at LO'"'^'''''' or in some hybrid 'nLO' approximation'''^' -'^^ 
in which the one loop anomalous dimensions are used together with one loop 
coefficient functions in some sensible scheme (one in which the first moments 
are given by (2.6)), in an effort to assess the effect of the direct coupling of Ag 
to gi . The calculation^^ of the two loop anomalous dimensions has now made 
complete NLO calculations^^ possible. 

At LO the effects of evolution are genc;rally small, and in particular the 
asymmetries are reasonably scale independent.^^ However at nLO the direct 
gluon coupling can have a dramatic effect, especially if Ag is large. In fig. la 
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we show the resuh of two nLO fits^"' to the proton data*'^" alone, one (the 
'minimal gluon') in which Agi{tQ) = 0, while Aqf{to) is fitted, the other 
(the 'maximal gluon') in which Ag^(to) = as {^^ expected from the Zweig 
rule) while Agi{to) = is fitted.'' The minimal gluon fit shows that there 
is considerable evolution from the E143 data up to the SMC data, but below 
the crossover (at x ~ 0.3) this is fairly uniform in x. The maximal gluon 
fit is much more dramatic: there is a second crossover at x ~ 0.03 while at 
yet smaller x g\{x,Q'^) becomes increasingly negative as increases. The 
evolution between the two crossovers is then correspondingly larger. This is 
particularly evident in the corresponding asymmetries (fig. 2a), which are fairly 
flat for the minimal gluon but rise quite steeply with for the maximal gluon. 

Although on the basis of these fits to the proton data alone it was not 
possible to distinguish between minimal and maximal gluon, the correspond- 
ing predictions fig lb for the deuteron structure function are more distinct, 
essentially because the deuteron is predominantly singlet, and a large direct 

^ In both cases Agj" is fixed by eqn.{2.3). Other details of the fits may be found in ref.l4. 
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Figure 3: The structure functions xg^(x, Q^) and xgf(x. Q^) plotted against x. The notation 
of the data points*"^^ and fitted curves'*^ is the same as in fig. 1. 

gluonic contribution to gi thus tends to make gf negative at small x. The 
data^'^^ seem to prefer a maximal gluon. The strong growth in the deuteron 
asymmetry which this generates is apparent in fig. 2b. 

This conclusion has recently been confirmed by a complete NLO cal- 
culation, with both proton and deuteron data*~^^ included in the fit.^^ 
The parton distributions turn out to be surprisingly well determined, with 
Aqf{to) = 0.5 ± 0.1, Agiita) = 1.5 ± 0.8. The Zweig rule expectation As ~ 
is thus confirmed experimentally, while the discrepancy in the Ellis- Jaffe sum 
rule^*^ is accounted for almost entirely by a direct polarized gluon contribution, 
just as was conjectured in ref.2. The behaviour of the parton distributions at 
large x is roughly consistent with quark counting rules. At small x the non- 
singlet distribution is singular (behaving as x~^, A = 0.7 ± 0.2) while the 
singlet quark and gluon distributions are generally either flat or valencelike, as 
expected from the theoretical considerations reviewed in the previous section. 

Using this fit (displayed in fig. 3) to determine the first moments F^, Tf 
and ao (as defined in eqn.(2.1)) we find^^ 

r? = 0.122 ± 0.013 (cxp.);;];°ii (th.), 

Tf = 0.025 ± 0.013 (exp.) tlZl (th-), (4-1) 
ao = 0.14 ± 0.10 (exp.)+°;i^ (th.), 

at = lOGeV^. The central values are lower than those given by the experi- 
mental collaborations because of the scale dependence both of the asymmetries 
in the measured region and of the small x extrapolations: the experimental 
uncertainty is larger because of the uncertainty in the size of the polarized 
gluon distribution which drives this evolution. There is also a theoretical error 
which is predominantly due to an estimate of NNLO corrections: these can be 
large because the two loop coupling of Agi to Fi is effectively NLO (compare 
(2.7) with (1.5) when N = 1). 
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5 Desiderata 

In conclusion, the effects of perturbative evolution have to be taken into 
account when extracting axial charges from polarization asymmetries, since 
they may be large due to the anomalous coupling to polarized gluons. Con- 
versely, the structure of the evolution seen in the combined proton and deuteron 
data sets indicates that the polarized gluon distribution may indeed be large, 
making a substantial contribution to the nucleon spin. In order to make this 
conclusion more definite, more data over a range of at moderate values of 
X from one experiment are needed. Polarization of the protons at HERA,^^ 
enabling a measurement of (and perhaps eventually A^\) would further pin 
down the small x contribution and the size of the gluon. Data for moderately 
small X (say 0.001 ^x^ 0.01) but over a wide range of would be most use- 
ful. Indeed the behaviour of polarized structure functions at small x promises 
to be a very rich subject, for both theorists and experimentalists alike. 
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